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We develop a model of string dynamics with back-reaction from both scaling and non-scaling loops 
taken into account. The evolution of a string network is described by the distribution functions of 
coherence segments and kinks. We derive two non-linear equations which govern the evolution of 
the two distributions and solve them analytically in the limit of late times. We also show that the 
correlation function is an exponential, and solve the dynamics for the corresponding spectrum of 
scaling loops. 
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I. INTRODUCTION 



o 
o 

Topological defects, such as cosmic strings, could have formed at the end of the brane inflation pj or during the 
symmetry breaking phase transitions [2J in early universe. It is also known that the cosmic string could lead to a 
number of observable phenomena: gravitational lensing 0, 0] , CMB non-Gaussianities [5|-[9| , gravitational waves [10| 
and ultra- high energy cosmic rays [111 ll 21 ] . However before one can effectively look for the observational signatures 
from cosmic strings it is absolutely necessary to know the statistical properties of the strings at late cosmological 
times. Unfortunately this problem proved to be highly non-trivial and after thirty years of numerical [13l - l21| and 
Q_i' analytical [22h28J studies the final word has not yet been said. 

Nambu-Goto evolution of strings is usually described by decomposition of the position three vector x(<7, t) into right 
r^ i ■ and left moving waves 

<3: X (a,t) = a(<J - t) + h{a+t) a) 

^ H ! . . . 

with condition |a'| = |b'| = 1, where prime denotes a derivative with respect to a. Without intersections the dynamics 
f***- ■ is completely linear in the flat space-time but some non-linearities already arise from cosmological stretching on the 

expanding backgrounds (3) ■ Another source of non-linearities is due to the production of scaling and non-scaling loops 
lf\ . [29| . In order to analytically describe the complicated dynamics of cosmic strings all of the non-linear effects must be 
f^ included in a self-consistent way. 

In the previous paper [29j we made a first step to develop a self-consistent model by including the back-reaction 

from only non-scaling loops whose sizes do not scale with time. The main assumptions of the model were: 

1) non-scaling loops are predominantly produced at the scales of the initial correlation length ~ l m in- 
"'~j ■ 2) large scale ~ t inter-commutations do not significantly affect the statistical properties on smaller scales ~ lmin- 

3) in the comoving coordinates l m in remains constant, when all other length scales grow linearly with time. 

4) the gravitational backreaction scale remains always smaller than the scale of initial conditions ~ I 



nun ■ 



The first three assumptions are motivated by many numerical simulations [16l 1 1 8142 11 ] as well as by analytical results 
[3, [25|,[26{, however, the latter assumption might eventually break down. The effect of the gravitational backreaction 
on the string network is unimportant for comparisons of our model with Nambu-Goto numerical simulation [la Il8l - l21j . 
but is very important phenomenologically and will be discussed extensively in a separate publication [301 ] . 

The main mechanism responsible for production of non-scaling loops is described in details in the Appendix. To 
summarize, if one chooses at random two opposite moving wiggles (one left-moving and one right-moving) of invariant 
length l m in, then the probability of a' and — b' curves to intersect (or to form a cusp) is given by 

p ~ (1 - C(l min ))/2 (2) 

where 

C(l) = (a'(O)a'(O) (3) 
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is a correlation function. In Ref. 29] it was shown that the statistical properties of the network are uniquely 
determined by what we called the number of directions: 

N=l/p~2/(l-C(l min )). (4) 

The non-scaling model predicts an appearance of two new length scales: the coherence length 

m = ^ ( 5 ) 



and the cross-correlation length 



X(t) = p (6) 



where c ~ k ~ 1, in addition to the well known inter-string distance and correlation length 

d(t) ~ C(*) ~ t (7) 

At the onset of evolution N ~ 10 and at late times N grows logarithmically 

N(t) ex log(t) (8) 

due to cosmological stretching and emission of small loops. 

In this paper we make one step further and develop a semi-scaling model with back-reaction from both scaling and 
non-scaling loops. In a fully scaling network the total string length decays as 

dm__ r L(t) r L(t) r L(t) 

j, J- small , J- larqe , -L friction , ■ \P J 

dt t t t 

where 

*■ small T 1 large \ J- friction — ^ V-^/ 

and the three terms describe the energy transfer into small loops, large loops and Hubble friction respectively. Ac- 
cording to [11]: 

t small ^ t large ^ -L \-^-^J 

and r friction *C 1. Thus, about half of the total energy goes into small loops and another half into large loops. 
Derivation of equations which describe both decay channels will be the main subject of this paper. In the second 
section we derive the evolution equation for the distribution of coherence segments and in the third section we derive 
the evolution equation for the angles between these segments. The main results are summarized in the conclusion. 

II. DISTRIBUTION OF SEGMENTS 

Consider the total number n(a, t) of coherence segments of length a at time t. Without loss of generality we assume 
that non-scaling wiggles have unit size l m in ~ 1 and from the first assumption (see above) most of the non-scaling 
loops are of a unit size as well. Then there are three mechanisms which can make the number of segments n(o~, t) 
change in time: 

1) partial decay of a segment into unit size loops, 

2) merger of a pair of segments into larger segments, 

3) complete decay of segments into larger loops. 

If we assume that on each unit time step only the leading unit size (left- and right- moving) wiggle is removed from 
randomly chosen coherence segments, then from ([5]), (j9|) and (fTTjl the probability for a given segment of size a to 
shrink to a segment of size ~ a — 1 is given by 

Pfrr 11^ Mil * ° <^\ 



where J2T=i n ( cr ' *) i s the total number of coherence segments and L(t)/t is the total number of unit wiggles to decay 
in unit size loops in unit time. 

We can now write down a finite difference equation with all three mechanisms taken into account: 

n(a,t + l)-n(a,t) = -n(a,t)- +fi((r + l,t)-+/jVj)(i,i)ti(M)-j)(ff- l-x,t)~ — (13) 

a;=l 

where 

, n(a,t) n(a,t)£(t) . . , .. ,. 

pa, t) = V '/ = V ; ^ W ex n (<r, i i 3 14 

is the probability for a random coherence segment to have length a, and \x is the probability of the next-to-nearest 
coherence segments to point in the same direction which is required for a merger to occur. Because any two next-to- 
nearest segments had passed through almost the same set of opposite moving segments they are very likely to point in 
the same direction: /i ~ 1. (For the dating toy model introduced in Ref. [29] the number was calculated numerically: 
/i ~ 0.58.) The first two terms on the LHS of (fTB")) describe a partial decay of segments, the third term describes a 
merger of two segments pointing in the same direction and the last term describes the decay into large scaling loops. 
By plugging (jj"4"]) into (flUf we get 

n(a, t+1)- n{a, t)~-^(n{a+ 1, t) - n{a, t) + -0-- J^ n(x, t)n(l, t)n(a - 1 - x,t)\ - ^^ (15) 

or in a differential form 

where C\ ~ C^ ~ 1. We can use (|14[) to write down an evolution equation for a probability distribution 



- 2 — V — — (- —r^p{l, t) / p(x, t)p(a -1-x, t)dx. (17) 



dt t N 2 da N 2 

Despite of its complexity the integro-diffcrantial equation (flT|) has a very simple solution in the limit 1 <C cr <C t: 

exp(-M^) 
P(v,t)<x t Clt (18) 

which could be checked by direct substitution. In the expanding universe there would be an additional effect which 
could lead to a merger of two coherence length segments and would effectively lead to a faster growth of the coherence 
length or equivalently larger C\. In fact this effect could be easily included in the non-linear term of (TT7|) with a 
proper redefinition of C2 , but the overall form of the solution (|18[) would remain unaffected. 

Although the exponential distribution of coherence segments (fT5j) does not directly imply that the correlations 
function is necessarily an exponential on all scales it certainly suggests that the latter is nearly linear on a wide 
range of scales between l m in and £, where only the linear term in the expansion of the exponential gives us the 
dominant contribution. This feature is a direct prediction of our model and is mainly due to constant mergers of the 
next-to-nearest coherence segments described by a non-linear term in (fP7|) . 

From (|18[) we can find the probability for a pair of unit size wiggles on a distance a to be belong to the same 

segment exp ( — 3 ^ ° 1 ~ exp ( — 7 ) which can be used to estimate the probability for any two wiggles on a distance 
a ^> £ to point in the same direction 

a\ /i / s / x \ ( a — £ — x 



P( ( r)^exp[--l+^y exp (-- I exp I ^ ) dx + ... (19) 

where the higher order terms have a similar structure. The resulting probability 

P(a) - exp U|J (1 +m|-..) - exp (- (1 ~ ^ ) (2D) 



but to get the correlation function one must sum over all N directions which would still give us an exponential under 
assumption of Gaussianity 

C(*)=acp(~), (21) 

where 

C ~ at (22) 

is the correlation length and a ~ 1. In what follows it will be convenient to consider a distribution of correlation 
segments a described by an exponential of (|21[) 

p(a,t) cxexp (-— J . (23) 

III. DISTRIBUTION OF KINKS 

Nearby coherence segments are separated by sharp kinks with angles 9 > \/yN. To understand the distribution 
of 9 we will derive an equation for the number density n(9,t)d6 in a sting network with total length in long strings 
Lit). There are three mechanism which can lead to changes in n(9,i): 

1) cosmological stretching, 

2) formation of loops, 

3) large scale inter-commutations. 
The first mechanism is described by 

d9 . . o, v 9 . 

— ex a(2(v 2 ) - 1)-, (24) 



where the scale factor grows as a(t) oc t a and 



2(v 2 )-l~-±, (25) 



according to [29(. This is due to the fact that the opposite moving wiggles cannot point in all directions, or otherwise 
the small loops would have been formed. Note, that the corresponding power-law decaying solution is valid only in 
the limit of large kinks 9 3> 1/yN. If the cosmological stretching would be the only mechanism, then we would write 
n{9 + d6, t + dt) = n(0, t) or in a differential form 

dn{9,t) a9 dn(9,t) 

-^r = C3 Ni^^- (26) 

By definition of coherence segments [29j , n(6, t) must vanish for small kinks 9 <C 1/y/N which can be imposed by an 
additional term 

dn(9,t) _ r< a9 dn(9,t) 1 n(9,t) 

~^T- C3 Nl^dB~~ CA WN~T- (27) 

where C3 and C4 are some constants of order one. 

The second mechanism is due to unequal chances for different kinks to survive the evolution. Clearly the segments 
connected by small kinks have much larger probability to remain untouched when large kinks can easily form loops 
by interacting with similar opposite moving kinks. The probability for a random kink to have size 9 is 

p(M = ^|Pocn(M)* 3 . (28) 

and the probability that in time £ it would meet an opposite moving kink of a similar size (on a logarithmic scale) is 
~ p(9, t)9. If these two kinks also intersect on a unit sphere (which happens with probability oc 9 2 ) then we are almost 
guaranteed to create a loop (see Appendix). We can now generalize (|27p to include the formation of intermediate 
loops 

d<9,t) a9 dn(9,t) 1 n(9,t) n(9,t) 2 a3 n(9,t) 

~^r~ = C3 ni^f- - C4 wn^~ + ° 5 ^w - ~r (29) 



where C5 ~ 1. The last terms comes from the assumption that the large loops can carry away random kinks without 
any preferences given to their angles 0. 

Large scale inter-commutations constantly inject sharp kinks into the network. If the segments which undergo such 
intersections are completely uncorrelatcd, then the angles of kinks will be distributed as oc sin(0). Clearly it is easier 
to create kinks with angle ~ 7r/2 when two random segments intersect with frequency estimated as oc L(t)/t 2 . The 
overall evolution equation with all three mechanisms combined is given by 

dn(8,t) a6 dn(8,t) 1 n(9,t) n(6,t)\ 3 L(t) n{0,t) 

-^tT = c "N-t^^' Ci ^N^~ C5 ^W "~ raW "T (30) 

where C§ is yet another constant of order unity. 
We can substitute J2H1) into d30l) : 



and look for a stationary solution p(9, t) = p{9): 

m.ojf- C lN m + C,rt W - C 8 1^ (32) 

where C3, C4, C5 and Cg are now some other constants ~ 1. 

The first term on the RHS of (|32|) introduces a sharp exponential cut-off for the smallest angles <C 1/yN and can 
be neglected for large angles S> 1/vJV. Moreover, in the limit of late times when N becomes significantly large the 
last term can also be dropped and we get a much simpler equation 

M$) _ n at p(0) , n w3 , 2 , m2 



,-ie 

For sufficiently large N the decaying solution goes as: 



CaN^y- + c^e^iey. (33) 



which can be combined with (118P to obtain 



p(0) « 1 (34) 



exp(-M^) 
p(g,(7,t)oc tg3 (35) 

in the limit t > 1, 1/-//V < < 1 and ? mm < cr < i, 

For large scales a ~ i and large kinks ~ 7r/2 the distributions p(er, t) and p(0, t) do not give an accurate description 
of evolution since the loops consisting of many coherence segments could also form. Instead it is convenient to study 
the evolution of the correlation segments p(a, t) and p(9, t), one of which was already derived (f2"3"|) without large-scale 
inter-commutations taken into account. In the limit of interest the large scale inter-commutations and formations of 
large loops are the main effects corresponding to the last two terms in (|32p . This equation can be rewritten for the 
angles between correlation segments 

d M = C 3 N^p{ef^C^ T ^. (36) 

Clearly, the distribution function for large angles does not depend significantly on and thus, we can estimate the 
spectrum of large loops (defined in [20]) by the distribution of correlation segments with = 7r/2: 

f(a,t)cxp(a,n/2,t)U-(j^ J (37) 

where 

d(t) - bt. (38) 

The second factor is due to the fact that large loops of size a could rejoin back to the network with probability 
proportional to their physical size ~ a 1 . 



The large kinks are produced by intercommutations of nearby strings with about one kink per time ~ d{t)/v per 
string length ~ £(£), where the average velocity v is given by (f^5|) . We can assume that the probability to create a 
large kink on a correlation segment p(a, n/2, t) is proportional to the area swept by the segment and to the distribution 
(f2"3"| generated by the dynamics on small scales. Then from ([57)1 and (|23[) the spectrum of loops is given by 




/(*,<)« j exp( - (-) |. C-iOi 



with a scaling peak at around a ~ ai ~ W. 



IV. CONCLUSIONS 



The main point of the paper was to derive the evolution equations (|17[) and pip with back- reaction from scaling and 
non-scaling loops taken into account. Remarkably the equations turned out to have very simple asymptotic solutions 
on a wide range of scales and sizes of kinks (|35[) . In addition, we obtained some preliminary results on the from of 
the correlation function (|2"TT) and on the spectrum of scaling loops (I3U1) . 

Throughout the paper we have explicitly assumed that the large scale inter-commutations and the small scale 
gravitational backreaction do not significantly affect the statistical properties on the scales of l m in- Although the 
assumption is very typical, it is not immediately clear why this should be the case in a non-linear system under 
investigation. Moreover it is already known that the opposite is not true and that the dynamics on the scales of l m in 
affects the evolution on all scal es f 291] . We are planning to study these and other related issues analytically and/or 
numerically in our future work [30(. 
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Appendix A: Small loops 

To understand the formation of small loops we can consider the following decomposition of the tangent vectors 



a'(cr) = a{a)x + P(a)y + ^l- a{a) 2 - /3{a) 2 z (Al) 

b'(ff) = -r(a)x + 5{<r)y -yjl- 7 (ct) 2 - 8(a) 2 z (A2) 

where x, y and z are the unit vectors. The above ansatz is quite general but nevertheless assumes that a' and b' 
have components along some direction (we call z axis) with opposite signs. This assumption is well justified for small 
wiggles in the vicinity of cusps (at a'(er) ~ — b'(cr) ~ z). In this limit the four functions (a, /3,7 and S) are much 
smaller than one, and by emitting higher order terms we get 



a'(o-) = a(a)x + /3(a)y + ( 1 - 



^i±mi , * (A3) 

b'(,l = ~(rn,- + ,»),/-( i-li?l±^L):, (A 4) 



A string self-intersects and a loop of size I can form only when 

a(<7 - <) + b(cr + t) _ a((j - £ + I) + b(a + 1 + 1) 
2 ~ 2 ' 

It is convenient to define three functions 



(A5) 



X(a,t) = °(*-*)+7(* + *) (A6) 

*>,*) - #*-*)+ '(*+*) (at) 



and write down the system of three integral equations 

f(T + l 



X(u,t)du = (A9) 

er + 2 

Y(u,t)du = (A10) 

-z 
Z(u,i)du = (All) 

with three unknown parameters a, I and t. It is clear that the solutions of the above system depends solely on the 
statistical properties of the three functions X, Y and Z which in turn are uniquely described by a, /3, 7 and 5. 

For the time being we can assume that the string is relatively straight on the scales ~ £ with only small wiggles 
with wavelength ~ l rnin and amplitude P(l m in)lmin- This means that all four functions a,/3,7 and S oscillate with 
wavelength ~ l m in and amplitude P(l m in) Cl. If we ignore the cross correlations between opposite moving waves, 
then both X and Y must also be oscillating functions in each of the variables (t or l m in) with similar wavelength and 
amplitude. However, the third function can be re-express as 

= (7(0- + t)- a(a - t)) Xja, t) + (5(a + t)- (3{o - t)) Y(a, t) _ 

which shows that the amplitude of oscillations is much smaller ~ P(l m in) 2 , but the wavelength is still the same ~ l m i n - 
In other words the variance on the scales l m in is given by 

^(X 2 ) - (Xp ~ v 7 ^ 2 ) - 00 2 ~ ^(U) and V(^ 2 > - (^> 2 ~ i'(fmm) 2 . (A13) 



8 

According to (|A12[) we are guaranteed to have Z(<r,t) — in the instant moments when both function X(a,t) and 
Y(a,t) vanish. In fact, one can always rotate the axis such that a(0) = ,8(0) = 7(0) = (5(0) = X(0, 0) = F(0,0) = 
Z(0, 0) = which corresponds to a mini-cusp at a = and t — 0. Since we are not interested in the zeros of X, Y and 
Z but in the zeros of their integrals we should integrate our three functions over a finite range I > as in (|A9|) , (|A10[) 
and (jAlip . In the limit of late times and small wiggles (l/£ <C P(l m in)) we can expect the mean of oscillations to 
change very slowly 

(X) ~ (Y) ~ (Z) ~ 0. (A14) 

which implies together with (|A13|) that the system of equations (|A9[) . (IA10[) and (|A11|) must have about one solution 
per region -l m in < t < l min , -l min < a < l min and -l mm < I < l min with I — t — a — l min . For l/£ < P{l m in) 
at least one mini-cusps occurs per length Z mi „ and therefore the entire overlapping region of the curves a' and —b' is 
likely to decay into small loops. 



